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Summary. In this study, we attempted to define differences in the hemagglutinin- 
esterase (HE) glycoprotein between 11 bovine coronaviruses (BCV) recent (post- 
1991) and past (pre-1991) isolates from neonatal calf diarrhoea (NCD) and winter 
dysentery (WD) syndromes as a basis for strain differentiation related to the clin- 
ical source of the isolates. The five WD-associated BCV isolates studied could 
be distinguished from past NCD-isolates by their hemagglutinating (HA) prop- 
erties at 4° and 37 °C, their receptor-destroying enzyme (RDE) activities with rat 
erythrocytes and lack of reactivity of these NCD isolates to four HA inhibiting 
(HAI) monoclonal antibodies (MAbs) directed against the HE glycoprotein of the 
reference WD-associated BCQ.2590 Quebec strain. Although minor or no dif- 
ferences could be demonstrated by comparing biological properties of the HE of 
WD.-isolates to those of recent NCD-isolates, past NCD isolates lacked reactivity 
with the WD HAI MAbs, whereas recent NCD isolates displayed two distinct re- 
activity patterns. Attempts to define sequence differences in the HE genes of the 
WD and NCD strains revealed high nucleotide (NT) identities with only scattered 
amino acid differences, seemingly unrelated to the clinical origin of the isolates 
or HAI MAb reactivities. 
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Bovine coronavirus (BCV) is an important enteropathogenic agent of the bovine 
species which replicates in the mature enterocytes of the small intestine and colon, 
causing severe diarrhoea in neonatal calves (NCD) [11, 16, 27] and chronic 
shedding in adult cattle (AD) [4, 12]. The virus has been reported to be also 
responsible for acute enteric infection in adult cattle during the winter season 
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[1, 8, 21]. The disease, usually referred to as winter dysentery (WD), is charac- 
terized by an acute onset of dark, liquid, hemorrhagic diarrhoea accompanied by 
decreased milk production and variable depression and anorexia in adult cattle. 
Other field isolates infect the upper respiratory tract from feedlot cattle leading 
to lesions of interstitial pneumonitis [3, 26]. 

BCV is a well characterized hemagglutinating coronavirus. The viral parti- 
cle is mostly spherical, enveloped, 120 to 160 nm in diameter, and displays a 
double fringe of surface projections [7, 8, 15]. The viral genome, comprised in 
an helicoidal nucleocapsid, consists of a large single-stranded RNA molecule 
with positive polarity, approximately 30 kb in length [24], and encodes the three 
major structural proteins of coronavirions: a phosphorylated nucleocapsid pro- 
tein (N) of 50 to 60 kDa, the integral matrix glycoprotein (M) of 20 to 38 kDa, 
and the spike (S) glycoprotein of about 200kDa which forms the long club- 
shaped surface projections [7, 24]. As other members of the hemagglutinating 
coronaviruses, including the hemagglutinating encephalomyelitis virus of swine 
(HEV), the human respiratory coronavirus (HCV-OC43), the murine hepatitis 
coronaviruses (MHVs) and the turkey enteric coronavirus (TCV), the BCV par- 
ticle also possesses an additional 140 kDa envelope-associated glycoprotein, the 
hemagglutinin-esterase (HE), which is made of two separate 65 kDa subunits 
linked together by disulfide bonds [7, 13, 15]. The HE glycoprotein is associated 
to the hemagglutinating and receptor destroying enzyme (RDE) activities [18, 19, 
25029). 

At present, only one serotype has been identified amongst BCV isolates, de- 
spite their involvement in different clinical diseases. Data obtained from recent 
serological studies, using polyclonal as well as anti-S or anti-HE monoclonal 
antibodies, suggest the existence of different subgroups of BCV, each containing 
strains from either NCD, WD or AD [11, 17, 18, 27]. On the genomic level, only 
minor variations in the S and HE genes have been demonstrated between virulent 
and avirulent or vaccine strains of BCV [3, 8, 30, 31]. Recently, the HE gene of 
one representative Quebec WD-associated strain (BCQ.2590) was sequenced and 
compared to that of the attenuated Mebus strain and two highly enteropahogenic 
NCD strains isolated during the same period [8]. Although 98% sequence ho- 
mologies were demonstrated with the Mebus strain, distinct aa substitutions were 
identified within the signal peptide and near the sequence of the putative esterase 
domain of the HE glycoprotein of the WD strain. Characterization of MAbs to the 
HE glycoprotein of the BCQ.2590 strain defined two distinct subgroups among 
BCYV isolates that have been associated either with typical outbreaks of WD or 
NCD in Quebec dairy herds from 1986 to 1996 [18] and variability was also 
demonstrated amongst NCD strains using anti-S MAbs [14, 17]. 

In the present study, we attempted to define biological, serological and ge- 
nomical differences in the HE between recent (post-1991) and past (pre-1991) 
BCYV isolates from NCD and WD syndromes as a basis for strain differentiation 
related to the clinical source of the isolates. 

The cell culture-adapted and attenuated Mebus strain (L9 vaccine strain) of 
BCV was obtained from the American Type Culture Collection (ATCC VR 874, 
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Rockville, MD). This prototype NCD strain was originally isolated from diarrhea 
fluid of a two-week-old calf [13] and attenuated following more than 30 successive 
passages in fetal bovine kidney (FBK) cells [16]. The BCV.Que strain, another 
attenuated NCD strain isolated in 1972 from diarrheic calves in Quebec, Canada, 
was propagated for at least 41 passages in Vero, followed by 10 to 15 passages 
in MDBK cells prior sequencing studies [5, 9, 19]. For the present study, the ref- 
erence Mebus strain, the Quebec reference WD-2590 strain [8], as well as other 
representative field isolates of BCV associated to either NCD or WD outbreaks 
obtained from non-vaccinated dairy herds from five geographical areas in South- 
ern Quebec [17, 18], were propagated for not more than 5 successive passages in 
human rectal tumor (HRT-18) cells in the presence of 10 U/ml bovine pancreatic 
trypsin, as previously described [7, 8]. Previous experimental infection studies 
demonstrated that WD and NCD field BCV strains inoculated orally to one- to 
two-week-old calves after only 5 passages in MDBK or HRT-18 cells were still 
highly virulent [20]. Infectivity titers of the cell culture-adapted BCV strains were 
determined by titration of clarified tissue culture medium using an end-point di- 
lution procedure and and calculaion of TCIDs9 s/ml, as previously described [5]. 
Hemagglutination activities of clarified infected culture supernatants (HA) and 
HAI tests with rat erythrocytes were also performed as previously described [7, 
18]. The extracellular virions of the various BCV isolates were purified from 
200 ml of clarified infected cell culture fluids by isopycnic ultracentrifugation 
on sucrose gradients [7, 8]. To determine the acetyl esterase (AE) activity of the 
various BCV isolates, aliquots (15 wl) of purified preparations, adjusted to con- 
tain approximately the same amounts of viral particles as determined by negative 
stained electron microscopy [6], were added to 1 ml of PBS containing 1 mM- p- 
nitropheny! acetate (PNPA). Hydrolysis of the substrate following an incubation 
of 5 min was monitored at a wavelength of 405 nm in a spectrometer (Spectronic 
2000, Milton Roy). Their receptor destroying enzyme (RDE) activities were de- 
termined, according to Storz et al. [25], and expressed as the ratio between HA 
titers determined after incubation at 1h at 4°C and HA titers determined after 
shifting to 37°C for 1 h, leading to activation of the enzyme activity and elution 
of the viral particles from the surface of erythrocytes [8, 25]. Four anti-HE MAbs 
(1D6-2, 1D6-3, 7B11-19, 9F2-1R) were obtained following fusion experiments 
with splenocytes from Balb/c mice which have been immunized with an enriched 
peplomeric fraction of the WD-associated BCQ.2590 strain [18]. The HE speci- 
ficities (65 and/or 140 kDa, corresponding to the monomeric and dimeric forms of 
the HE glycoproteins) of the MAbs were previously determined by Western blot 
and immunoprecipitation [18]. No reactivity was observed toward the S (100 or 
200 kDa) glycoprotein. The MAbs were either of the IgG1 or IgG2b subisotype. 

The infectivity titers of the various field BCV isolates ranged between 10°°- 
10’ TCIDso s/ml after 3 to 5 successive passages with HA titers at 4°C varying 
between 1:320 to 1:2560 (Table 1). In general, the WD isolates displayed higher 
AE activity than NCD isolates, with the exception of NCD isolates BCQ.701 
and BCQ.1523. The five WD isolates tested had RDE activity titers that varied 
from 2 to 16 after 5 passages in HRT-18 cells. However, in the cases of WD 
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strains BCQ.2590, BCQ.2508 and BCQ.7373, higher RDE activity titers (up to 
256) were detemined at lower cell culture passages (first two passages) and thus 
were reported as ranges in Table 1. All NCD-associated BCV isolates including 
the Mebus strain, with the exception of BCQ.701 (RDE titer of 8), appeared to 
possess a minimal RDE activity since no reduction in their HA titres was noticed 
when shifting to an incubation at 37 °C. 

Data of cross-reactivity of MAbs directed against the HE of the WD-associated 
BCQ.2590 strain to the other BCV isolates analysed are also depicted in Ta- 
ble 1. With those MAbs, major differences could be demonstrated between the 
HAI reactivity patterns of the WD-associated BCV isolates and those obtained 
with the NCD-associated isolates. Indeed, the four anti-HE MAbs displayed 
similar HAI antibody titers than that obtained for the homologous WD.2590 
strain with the four other WD-associated isolates tested, except that MAb 9F2- 
1R failed to react to isolate BCQ.2439. However, none of these anti-HE MAbs 
showed reactivity toward the reference Mebus strain, as it was also the case with 
other recent (pre-1991) NCD isolates BCQ.571 and BCQ.3, but 2 distinct pat- 
terns of HAI MAb reactivity was observed with NCD isolates collected after 
1991. Indeed, NCD isolates BCQ.3708, BCQ.376 and BCQ.701 behaved like the 
older NCD isolates, whereas BCQ.1523 had a reactivity profile identical to WD 
isolates. 

For the preparation of genomic RNA, aliquots (50 wl) of purified virus or in- 
fected cell-lysates were supplemented with 1 wl of RNAguard (Pharmacia) and 
RNA extraction was performed by the one-step guanidinium isothiocyanate-acid 
phenol method [2]. The HE coding regions were amplified by RT-PCR as already 
described [8] using the following primers: oligo HE.B (5’-TTAATCTCCTATAA- 
CAGC-3’) an antisense primer representing the sequence at position 40 to 58 of 
the S protein gene [19] and, oligo HE.A (5’-GGTTTTATGAATCTCCAGTTGA- 
3’), a sense primer corresponding to the sequence — 43 to — 25 nt upstream of 
the 5’ end of the HE gene [19]. Following purification using the GeneClean 
II nucleic acid purification kit (BIO 101, LaJolla, CA), the RT-PCR amplified 
products with A overhangs were ligated into a TA cloning vector (pCRII vector, 
Invitrogen Co., San Diego, CA), providing single 5’ T overhangs at the inser- 
tion sites [8]. Sequencing of cDNA clones was performed on both strands by the 
dideoxynucleotide chain-termination method [22], using T7 DNA polymerase 
(Pharmacia Biotech Inc.) in an Automated Laser Fluorescent DNA sequencer 
(Pharmacia LKB). Sequences and additional primers used for sequencing HE 
gene were: oligo $418 (5’-AAGAATATGGCTGTGTATCG-3’), oligo S881 (5’- 
TTGATTCACGGTGGAACAATG-3’), oligo $1083 (5’-TGATAATGTTAGCA- 
GTGTCTGG-3’), and the four corresponding reverse complementary primers. 
To assess the error rate of the reverse transcriptase and Jaq polymerase, clones 
from 3 different RT-PCR events were sequenced. Errors generated by T7 poly- 
merase were avoid by sequencing randomly selected clones in both directions. 
Subsequently, the nt and deduced aa sequences were computer analysed with the 
aid of the Mac Vector 3.5 (International Biotechnologies) and Gene Works 2.4 
(IntelliGenetics inc., Mountain View, CA) programs. 
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The HE genes of four additional Quebec WD-associated BCV field isolates, 
and four antigenically distinct NCD-associated BCV field isolates were cloned 
and sequenced. The nucleotide sequence accession numbers (EMBL/GenBank/ 
DDBJ) that have been assigned to the HE genes of these Quebec enteropathogenic 
BCV isolates are AF230523 for BCQ.2508; AF230524 for BCQ.2442; AF230525 
for BCQ.2439; AF239306 for BCQ.7373; AF239307 for BCQ.1523; AF230526 
for BCQ.376; AF230527 for BCQ.701 and AF230528 for BCQ.3708. Compari- 
son study was made with sequences of the HE genes of the virulent and attenuated 
(L9) reference Mebus strain (GenBank accession numbers M30612 and M76372) 
and the highly virulent NCD strain LY 138 (accession number M76374) [30], both 
isolated in the early 1970’s, the BCV.Que strain associated to an acute outbreak of 
NCD in a Quebec dairy farm in 1972 [19], in addition to those of the previously 
sequenced Quebec WD-associated strain BCQ.2590 isolated in 1991 (accession 
number L38962), and Quebec NCD-associated strains BCQ.571 (accession num- 
ber V06093) and BCQ.3 (accession number L38963) isolated in 1985 and 1989, 
respectively [8]. In all cases, the HE gene was 1,272 nt long encoding a predicted 
protein of 424 aa residues (data not shown). Frameshift, deletion or insertion, 
and nonsense mutations were not observed. Nucleotide substitutions affected 45 
positions of the nt sequence among the BCV isolates studied, of which 30 (66.7%) 
involved purine substitutions (C—T; T-C), occuring mainly at wooble bases, thus 
leading to silent aa mutations (Table 2). The HE gene sequences of the four ad- 
ditional WD-isolates contained approximately 22 nt substitutions in comparison 
with the reference BCQ.2590 strain, a comparable number of nt substitutions be- 
ing identified in the cases of the four antigenically-distinct NCD strains analysed 
when compared with both the Mebus and BCQ.2590 strains. The percentages of 
nt identity between BCV isolates associated with both enteric syndromes var- 
ied between 97% to 100%. The four additional WD-associated BCV isolates 
possessed 99% nt identity with the reference BCQ.2590 strain and the highly vir- 
ulent NCD-LY 138 strain. The WD-associated isolates BCQ.2442 and BCQ.2508 
showed only two nt substitutions between them, and both showed only three nt 
substitutions when compared to the HE gene of the recent NCD-associated isolate 
BCQ.1523 which displayed similar reactivity profile to the four anti-HE MAbs. 

Data from the comparison of the deduced aa sequences of the eight BCV field 
isolates studied, with those of the previously sequenced Quebec strains isolated 
from 1977 to 1986 (BCV.Que, BCQ.3 and BCV.571), and those of the WD- 
associated strain BCQ.2590 and NCD-associated reference strains Mebus and LY- 
138 are summarized in Table 3. A total of only 17 aa substitutions were identified 
amongst the field isolates analysed in comparison to the deduced HE protein 
of the Mebus strain, which remained highly conserved following attenuation by 
numerous passages in cell cultures (L9 strain) [19, 31]. 

Analysis of the predicted aa sequences of the HE of BCV isolates studied also 
revealed that the 9 putative N-linked glycosylation sites (NXS/T) previously iden- 
tified for the NCD-associated Mebus and LY-138 strains, and the WD-associated 
BCQ.2590 strain, located at aa positions 54, 88, 153, 246, 301, 316, 358 and 
417 [8, 30], were conserved in BCV field isolates associated to both enteric 


Table 2. Nucleotide substitutions in the HE gene of Quebec NCD- and WD-associated field BCV isolates 


as compared to the Mebus, BCV. Que and LY-138 reference strains 


BCYV field isolates or reference strains 


WD WD WD NCD 


Nt 


NCD LY BCV Meb 


NCD NCD NCD 
701 571 


NCD 


location WD WD 


L9 


138 Que 


3708 3 


2590 2508 2442 2439 7373 1523 376 


RRR CRUST ORRUDTURUURURHORURHRERHOUUSTELOUSE HR ROUUUUUER 


* * * * 


* * x * 


* 


* 


x» Ox x x 

Ox « Ox « # « » 

Ox Bx « Ox « # 

Ox « » « Ox x « 

Ox « « « Ox x x 

Ox « *« BOx « « 
xe « Lx 


* 


* 


* * * * * * * * * * * * * * 


ORR« UUSER« Ux UV« 


ORR« UUSER« Ux UV« 


xx» Bx OULTRUVUOVUD: 


* * * * * 


x Ex * * 


BCV Bovine coronavirus; Meb NCD reference Mebus strain; WD winter dysentery; NCD neonatal calf 


diarrhea Meb.L9 and BCV.Que are attenuated BCV strains; all others isolates, including LY-138, are virulent 


strains. Nt location given according to Zhang et al. [27] 
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Table 3. Amino acid sequence substitutions in the HE glycoprotein of Quebec NCD- 
and WD-associated field BCV isolates 


AA Quebec BCV field isolates and diseases BCV reference strains 


location WD WD WD WD WD NCD NCD NCD NCD NCD NCD LY BCV Meb 
2590 2508 2442 2439 7373 1523 376 701 3708 3 571 138 Que L9 


5 P P P P P P * * P P * Px L 

8 A * * * * * * * * * * * * Vv 

11 * C C * * C * * * * * * * S 
49 * T T * T T T T T T T * * N 
53 P * * * * * * * * * * * * R 
66 * * * * * * G G G * * * * D 
103 * * * I * * * * * * L * * Vv 
114 * * * * * * * * I * * * * T 
182 * * * * R * * * * * * * * G 
245 * * * * L * * * * * * * * F 
282 * * * * G * * * * * * * * D 
344 A * * * * * * * * * * * * D 
350 R * * * * * * * * * * * Ok Q 
367 P P P P P P P P P * * Px S 
379 * * * * * * * * * * * Nx T 
391 * * * * I * * * * I * * * P 
392 * I I I I I I I I I I * ox L 


BCV Bovine coronavirus; Meb reference Mebus strain; WD winter dysentery; NCD neonatal calf diarrhea 
Strains Meb-L9 and BCV. Que are attenuated strains; all other field WD and NCD strains, including LY-138, 
are virulent. There were a total of 424 amino acid residues for each HE gene sequenced. The first column 
refers to the Quebec WD reference strain BCQ.2590, as it differs from the reference BCV Mebus strain in the 
final column. The asterisks denote consensus with the BCV Mebus strain. A total of only 17 aa substitutions 
were identified amongst the field isolates analysed in comparison to the Mebus strain. 

The AA locations on the HE glycoprotein were given according to Zhang et al. [27] 


syndromes, as well as the putative esterase active domain FGDS (aa positions 36 
to 39). However, similarly to what was previously described for the highly virulent 
NCD-strain LY-138 and WD-strain BCQ.2590 [8, 31], all four additional Que- 
bec WD-associated isolates, as well as the NCD-associated isolates BCQ.1523, 
BCQ.3 and BCQ.3708 contained a proline substitution at aa position 5 of the 
predicted signal peptide, which was not the case for the other NCD-associated 
isolates studied. It has been previously postulated that this substitution may al- 
ter the maturation and intracellular transport of the HE glycoprotein, which still 
remained to be elucidated [30]. It also appeared that all WD-associated BCV 
isolates possessed the proline substitution observed for BCQ.2590 and LY-138 
virulent strains at aa position 367, as it was also the case for the Quebec NCD- 
associated field isolates collected after 1990. Further analysis revealed high NT 
identities with only a few scattered aa differences, seemingly unrelated to the 
clinical origin of the isolates or the HAI MAb reactivities. 
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As previously reported for the reference Quebec WD-associated BCQ.2590 
strain [8], data obtained in the present study showed that the four additional 
WD.-associated field isolates analysed were different from the prototype BCV. 
Mebus strain in respect to clinical symptoms, HA patterns at 4°C and 37°C, 
and higher RDE activities with rat erythrocytes. However, similar HA and RDE 
properties were observed for two of the recent NCD-associated isolates studied 
(BCQ.1523 and BCQ.701), meaning that these characteristics are not unique to 
WD.-associated isolates, but at least BCV isolates associated to the attenuated 
Mebus (L9 vaccine strain) [16, 25] and BCV.Que [5, 19] strains can be distin- 
guished by their capacity to hemagglutinate at 37 °C and their lower RDE activities 
when tested with rat erythrocytes. At the molecular level, two unique Pro substitu- 
tions were previously identified for the HE of the reference WD-strain BCQ.2590, 
one in the signal peptide (aa position 5) and one near the putative esterase do- 
main FGDS (aa position 53) [8]. Interestingly, the FGDS sequence was found to 
be well preserved in all BCV isolates studied, as it was also previously demon- 
strated by others by comparing virulent and avirulent BCV strains, MHV-A59, 
MHV-JHM, and influenza C viruses [8, 19, 28]. The data obtained thus indicate 
that the esterase active site is essential for the structure and function of the HE gly- 
coprotein. However, as it was previously hypothesized, this domain is probably 
not a determinant for BCV virulence, since it is conserved in both virulent and 
avirulent strains, as well as in strains associated to various clinical symptoms [3, 
30]. Whether the Pro substitutions at aa positions 5 and 367 are involved in viral 
infectivity, as previously suggested [3, 8, 31], is still to be demonstrated. Inter- 
estingly, both substitutions were present in the HE protein of all WD-associated 
BCV isolates tested, and both NCD-associated isolates BCQ.1523 and BCQ.701 
showing higher AE activities. 

Although there exist BCV strains that can be distinguished by their clinical 
syndromes, the existence of distinct BCV serotypes has yet to be elucidated. 
Previous investigators detected only minor strain variations amongst cell culture- 
adapted BCV strains isolated from NCD cases in Europe, and from enteric and 
respiratory cases in the United States by both IIF and HAI tests [11, 27]. More 
recently, anti-S MAbs were found to distinguish between vaccine (avirulent) and 
wild-type (virulent) BCV strains in mildly denaturating Western immunoblotting 
and neutralization assays [14]. We also previously reported antigenic plurality 
amongst BCV isolates from NCD cases in Quebec through the partial reaction in 
an indirect ELISA of MAbs against the S glycoprotein [17]. In the present study, 
the four anti- HE MAbs which inhibited HA activity of the homologous BCV.2590 
strain did not react to the reference Mebus strain, while displaying strong HAI 
titers to the other WD isolates tested. In disagreement with our previous findings 
[18], the reactivities of MAbs 1D6-2, 1D6-3, 9F2-1R and 7B11-19 were not 
restricted to WD strains as similar HAI titers were obtained with one of the recent 
NCD- isolates (BCQ.1523) tested. Although data demonstrated that variations 
affecting antigenic determinants of the HE glycoprotein occurred among BCV 
strains isolated in Quebec within the last ten years, this diversity appears to be 
unrelated to the clinical symptoms, which is in agreement with the recent findings 
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by Tsunemitsu and Saif [27] in the case of BCV strains isolated in the United 
States. 

Although there seems to have been an alteration in antigenic determinants of 
the HE glycoprotein of Quebec isolates after 1991, the sequence analyses did not 
show many distinct variations that are common to all field isolates, which would 
distinguish specific epitopes targeted by each of the MAbs. One may postulate 
that the MAbs studied are directed toward conformationally dependent rather than 
linear epitopes, but a conformation change of the protein and not lack of epitopes 
could also explain the serological data obtained. In a previous study [18], it has 
been demonstrated that two of the anti- HE MAbs tested (MAbs 9F2-1R and 1D6- 
3) failed to react by Western blot toward the HE protein of the homologous strain, 
suggesting that these MAbs are more probably directed toward conformationally 
dependent epitopes. However, conformation changes of the protein could have 
resulted from aa variations, such as the substitution of Pro residues in certain 
regions which could entail the alteration of the three-dimensional structure of the 
protein, thus serving to either expose or mask certain epitopes present in some 
strains and not in others. Similar conclusions may be drawn about the higher RDE 
activities of various strains, as a change in conformation of the protein could result 
in an enhanced enzymatic activity through a more prominently exposed region. 
Sequencing analysis of more strains that have been isolated prior 1991 are needed 
to support this statement. 

In previous studies, HAI was also the only serologic test hat allowed differ- 
entiation between BCV strains associated to both enteric syndromes [1, 8, 11, 
27]. Since only few nt substitutions could be revealed in the HE genes of the WD 
and NCD BCYV strains, such serological differentiation may rather involved the 
S glycoprotein which has been also associated with HA activity [23]. Retrospec- 
tive studies of the humoral immune response to individually expressed recombi- 
nant HE and S glycoproteins are required to determine which of the two major 
envelope-associated glycoproteins is mainly involved in the antigenic variability 
of recent and older NCD- and WD-associated Quebec BCV field isolates. 
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